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Abstract. We discuss various aspects of vector meson production, first analysing the interplay between 
perturbative and nonperturbative aspects of the QCD calculation. Using a general method adapted to 
incorporate both perturbative and nonpertubative aspects, we show that nonperturbative effects are im- 
portant for all experimentally available values of the photon virtuality . We compare the huge amount 
of experimental information now available with our theoretical results obtained using a specific nonpertur- 
bative model without free parameters, showing that quite simple features are able to explain the data. 
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1 Introduction 



Electroproduction of vector mesons provides an interest- 
ing laboratory for studying the interplay between per- 
turbative and nonperturbative QCD. Most emphasis has 
been put on the perturbative side of the calculations [TJ 
EUSEIISIEIIIIIHIEIIIO]. There the production process is con- 
sidered to be mediated by two gluon exchange and the 
coupling of the gluons to the hadron, generally a proton, 
is taken from a phenomenologically determined gluon dis- 
tribution as obtained from deep inelastic scattering (DIS). 
The coupling of the exchanged photons to the produced 
quarks is treated perturbatively. This is justified if there 
is a truly hard scale and the produced quarks stay close 
together. 

The emphasis put on perturbation theory is understand- 
able given its justification in terms of basic elements of 
QCD. Nevertheless we think it is useful and even necessary 
to scrutinize also the other side of the medal, namely the 
nonperturbative aspects, and investigate its consequences. 
To examine the roles and magnitudes of different effects 
we have chosen an approach which starts from general 
expressions and investigate the limit where perturbation 
theory holds. We derive and discuss particularly the devia- 
tions from perturbative QCD induced by nonperturbative 
effects. 



— Even if the produced vector mesons are heavy, the ef- 
fect of binding by a confining potential is not negligible 
at presently accessible values of the photon virtuality 

. The binding effect influences very strongly the 
dependence of the production amplitude. 

— Even for high values of the production of trans- 
versely polarized vector mesons receives important con- 
tributions from regions where the two produced quarks 
are widely separated. 

— In spite of a clear transition from the perturbative 
to the nonperturbative regime, there are nevertheless 
striking systematic features in the production of light 
and heavy vector mesons, whose understanding requires 
the incorporation of nonperturbative methods. 

— The gluon distribution of the proton as obtained from 
DIS allows only to calculate production amplitudes 
with zero momentum transfer. This requires quite es- 
sential extrapolations in the analysis of the experimen- 
tal data. Since the approach discussed here is based 
on a space-time picture it takes skewing effects auto- 
matically into account, and therefore the dependence 
of the production amplitude on (moderate) values of 

the momentum transfer from the virtual photon to 
the proton, can be calculated. 

— In QCD the production of vector mesons is closely re- 
lated to purely hadronic scattering processes without 
requiring the use of vector dominance models. In order 
to obtain a unified picture, a handling of nonperturba- 
tive effects is clearly necessary. 



There are several important reasons for this approach. 
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— Closely related to the previous item is the relation 
between QCD and Regge theory. In order to inves- 
tigate possible bridges between the two approaches 
again nonperturbative methods must be used. 



The calculation of photo and electroproduction presented 
in this paper uses a general method for high energy scat- 
tering based on the functional integral approach to QCD 
and on the WKB method, which is capable of incorpo- 
rating both the perturbative and nonperturbative aspects 
|lipi2j. The nonperturbative input is given by a special 
model of nonperturbative QCD, called stochastic vacuum 
model [TTiTi] , that has been successfully applied in many 
fields, from hadron spectroscopy to high energy scatter- 
ing. The very satisfactory results of the model in purely 
nonperturbative regions gives a strong weight to our de- 
termination of nonperturbative correction terms near the 
perturbative regime. The energy dependence is based on 
the two-pomeron model of Donnachie and LandshofF [15j . 

In previous papers we have investigated photoproduction 
[TB] and electroproduction [17] of 3/ip vector mesons. We 
have also used the same framework to investigate some 
general features of photo and electroproduction of all S- 
wave vector mesons |18j which arise from the structure of 
the overlaps of photon and vector meson wave functions. 

In the recent years many more data have been obtained in 
HERA experiments, with higher accuracy and statistics, 
and their comparison with theoretical calculations provide 
opportunities to understand and describe important gen- 
eral features of the dynamics governing these processes. 
In particular, we may learn in what extent the exper- 
imentally observed features are contained in the global 
nonperturbative aspects of the systems in the final and 
initial states, such as their wave functions and the long 
range correlation properties of the intervening forces. In 
this paper we describe very successfully most of these re- 
cent data, using the same framework that has been tested 
in several other cases, without the introduction of any free 
parameter. 

Our paper is organized as follows. 

In Sec. 2 we describe the methods used in the theoreti- 
cal calculations of photo and electroproduction of vector 
mesons. Since this has already been done on several oc- 
casions, we only give a short and schematic description, 
and then focus the attention on the comparison with the 
usual perturbative treatment. 

In Sec. 3 we present the results of our calculations of photo 
and electroproduction of all IS- wave vector mesons and 
compare them with the experimental data. We show that 
the predictions for all observables, the absolute values as 
well as the dependence on the photon virtuality Q^, the 
momentum transfer t and the energy W are very satisfac- 
tory. 

on our results and present them in context. 

The appendix gives formulas and details of the theoretical 
calculations and presents some of their general properties. 




transverse collision plane 



Fig. 1. Loop-loop scattering. 

2 General formulae 

and the strictly perturbative limit 

We start from a general approach to scattering based on 
functional integrals and the WKB approximation pTlfT^ , 
which has been adapted to hadron hadron scattering [ini 
|20] and to photo and electroproduction of vector mesons 
hadrons pTl fTOlfTT] . The basic features can be seen in Fig.[l] 
that represents the loop-loop scattering amplitude in real 
space time. The space time trajectory of the photon is rep- 
resented by a quark-antiquark loop, that of the proton by 
a quark-diquark loop. The transition to observable elec- 
troproduction amplitudes of hadrons is achieved through 
a superposition of the loop-loop amplitudes with the light 
cone wave functions of the hadrons and the photon used 
as weights. This leads to an electroproduction amplitude 
of the form 

T^^p^VpAW, t; Q2) = {-2iW^) (1) 
j SRidzii!v\{zuRi)*il^-y'x{zi,Ri,Q'^)J{q,W,zi,Ri) , 

with 

J(q,I¥,zi,i?i)= (2) 
j d''R2<fbe-"i'' |Vp(i?2)|' S{b,W,z^,R^,Z2^l/2,R2) . 

Here S{b, W, zi,Ri, 1/2, R2) is the scattering amplitude of 
two dipoles with separation vectors Ri , R2 , colliding with 
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impact parameter vector b; q is the momentum transfer 
vector of the reaction, with 



In these expressions = —Pj* is the photon virtuahty, 
W is the center of mass energy of the 7*-proton system, 
and t — {py — p^*Y is the invariant momentum transfer 
from the virtual photon to the produced vector meson; 
zi is the longitudinal momentum fraction of the quark 
in the virtual photon and in the vector meson, and we 
call zi = 1 — zi; ip^ZjR) represent the light cone wave 
functions. 



The differential cross section is given by 

da _ 1 
d\t\ ~ 167rH/4 I 



in 



(4) 



For the special case of forward production (q = 0) this ap- 
proach reduces to the dipole model for electroproduction 

m- 

Details of the evaluation of the loop-loop amplitude in 
the stochastic vacuum model can be found in previous 
publications [TOll^ . 



2.1 The perturbative limits 
and nonperturbative corrections 

For the photon wave functions ip^t\{zi, Ri,Q'^) we use 
the well known lowest order expressions, keeping the same 
notation used before p!Blll7j . 

If we concentrate on high values of virtuality the two 
valence quarks stay close together and we may use the 
dipole cross section of a small object of size Ri , that is 



J{q = 0, W, zi,Ri) = C{W, Q2) Rl 



(5) 



If we furthermore ignore the transverse extension of the 
vector mesons, replacing 

V'yA(zi, ^pv\{zl,0), (6) 

we can perform the i?i integration in Eq. ^ explicitly. 

Later we introduce specific ansatze for the vector-meson 
wave functions ipv\{zi, Ri) and use the results of a specific 
model, the stochastic vacuum model, for the evaluation 
of the loop-loop scattering amplitude J{q, W, zi, i?i), but 
for the moment we stick to the general formulae and first 
study the corrections to the strictly perturbative limit. 

We first consider the simpler case of longitudinally po- 
larized photons, where we obtain with the replacement 
(|6| (for later discussion we keep the dependence on the 
quark mass to/ in our expressions): 



npert 
VSQ! 



(7) 



Q I dzi 



Stt zjzf 



{zizii 



The wave function at the origin ^pvo{zi,0) is related to 
the coupling of the vector meson to the electromagnetic 
current fy by 



V3 



dzi 162;i2;i'(/)vo(2:i, 0) 



(8) 



Here e/ is the quark charge, and ey is the effective charge 
in the meson, that is ey = ly/2 for the p and ey — l/3\/2 
for the Lo meson , while for the (j), ■0 and T mesons we have 
By — Bf. 

We introduce 



»i 



dzi 







izfzf tljyo{zi,0) 

(ziZlQ^ -I- TO^)2 ; 



(9) 



dzi zizi^yo{zi,0) 



and then write the perturbative expression as 



(10) 



TT ey 



Q 



(g2/4 + m2)2 



Noting that the longitudinal wave function is suppressed 
at the end points zi — 0, zi = 1, we see that in the limit 
Q2 _> oo the mass to/ can be neglected against ziZiQ^ 
and the expression rj^ becomes independent of 



/q dzi ^pyo{zi,0) 
4 Jq dzi ziZiipvo{zi,0) 



(11) 



This is the correction factor due to the longitudinal ex- 
tension of the meson as obtained in [3 (up to a factor 2 
due to different definitions). However, for practical pur- 
pose we keep using the expression ( 10 1 with finite quark 



masses. If there were no binding effects, zi would be 1/2 
and therefore Ty^ would approach 1 for going to infinity. 



Comparing Eq. (Ill with the perturbative expression for 
the longitudinal electroproduction of vector mesons [3] we 
obtain 



CiW, 



(12) 



where x is the Bjorken variable x = (Q^ + M^)/W'^. This 
yields the well known relation between the dipole cross 
section and the gluon density |23j . After introducing our 
spe cific model for the meson wave functions in subsection 
2.3 we will present numerical values for ?7l((5^), which rep- 
resents the longitudinal momentum correction to the pure 
perturbative calculation of the amplitude for longitudinal 
photons. 

Nonperturbative effects also lead to a finite extension of 
the vector meson and to a deviation of the simple quadratic 
behaviour of the "dipole cross section" J{q = 0,W, zi, Ri) 
in Eq. ^ . We present numerical values for the rather large 
corrections due to these effects in subsection 12.31 



As it is well known, the treatment of the transverse po- 
larisation is more delicate for at least two reasons. 
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There is no strong suppression of the photon wave 
function at the end points Zi = 0, zi = I and there- 
fore the effective scale, namely zizi might be quite 
low even for highly virtual photons. This makes among 
other things a systematic analysis impossible since 
the factor corresponding to Eq. (Ill diverges due to the 



singularities at the end points zi = 0, 1. 
— The meson wave function is supposed to be more com- 
plicated and relativistic contributions are likely to be 
important even for rather heavy mesons. 

For simplicity we assume that the wave function of the 
vector meson has the same tensor structure as the trans- 
verse photon. Then the overlap function ^ brings with 
the replacement ^ the form 



rppert 
^"p — >Vp.\—l 



-2iW^C{W,Q^)ef 



'6a 



27r 

;2^„/,(l)/ 



(13) 



/o 
where 



\ 4m^ ^yi(zi, 0) + 16 uj-^izj + zl)^);[{z^,Q) 



(zizi( 



-I 



d 



4'vi (2:1,0) 



2cj2 \diRl) 

The relation to the decay constant fv is here given by 
/6 1 



(14) 



fv 



I 



ev 



Mv ^/47r 



(15) 



1 m) i;vi{zuO) + 2uj''{zl + zl)^^^l{zi,Q) 
dz\ 



This leads to a relativistic correction factor 77T which is of 
order 0{uj^ /mi), given by 



2^2 



^t(Q') = (Q74 + ™^) 



(16) 



1 4m^ V^vi(zi, 0) + 16 u;2(zf + zl)^^y\{zx,Q) 







(ziZiQ2+to2)2 



1 4m2^v-i(zi, 0) + 8u;2 (z? + z\) ^'[.'^(^i, 0) 

(1Z\ ^ 

4zizi 



The final result for the for transverse polarisation is writ- 
ten analogous to Eq. (Ill, that is 



mpert 
7*p — A— 1 



^t(Q') 



Q 



C{W, Q2) 



(g2/4+m2)2 



(17) 



Only in the absence of binding effects it would he lu — 
0, zi — i, and we would then have ririQ^) — > 1 for 
00, as it is in the case of longitudinal polariza- 
tion. Numerical values for the correction coefficients rjT, Vl 



evaluated in a specific model are given in subsection 2.3 



Strong deviations of 77^ and rjx from unity indicate im- 
portant deviations of the amplitudes from the pure per- 
turbative calculation, due to longitudinal momentum in 
the photon-meson overlap and to relativistic effects. 

Important additional corrections arise from the finite ex- 
tension of the vector mesons, and they are also discussed 
in subsection 12.31 



2.2 The specific non-perturbative model 

In order to obtain numerical information on the correc- 
tions discussed in the two previous subsections, we need 
to use models both for the wave functions and for the 
interaction of gluons and hadrons. 

The photon enters through its usual perturbative wave 
function. In model calculations [H] it has been shown that 
the perturbative wave function can be used also for small 



values of Q 



f 



if an appropriate constituent mass is 



introduced. For light mesons, the masses have been de- 
termined by comparison with the phenomenological two 
point function for the vector current. The approach has 
been successfully applied in the theoretical calculation 
of structure functions, Compton amplitudes and photon- 
photon scattering. 

The wave function has been adapted from photons to vec- 
tor mesons including relativistic corrections motivated by 
the structure of the vector current [r7 l [2l ] [24 j [25] . As in 
previous papers |161ll7j , we use two forms of meson wave 
functions: the Bauer-Stech-Wirbel (BSW) [55], and the 
Brodsky and Lepage (BL) [27] forms, which are detailed 
in the Appendix. We use this type of wave function for all 
vector mesons, with quark masses determined from a best 
fit to the vector current (they are m„ = rjid = 0.2, rrig = 
0.3, TOc = 1.25, rrib = 4.2 GeV) . 

Our nonperturbative treatment of the high energy pro- 
cess is based on functional integration PTIIT?] and on the 
stochastic vacuum model [13,14]. The method has been 
described in several occasions [20 , 28','25] , and we only quote 
here a few of its characteristic features. The model is 
based on the assumption that nonperturbative QCD can 
be approximated by a Gaussian process in the coulour 
field strengths; the gluon field correlator is therefore the 
quantity determining the full dynamics. Its parameters 
are taken from lattice calculations [3D] . The model yields 
confinement in non-Abelian gauge theories and leads to 
realistic quark-antiquark potentials for heavy quarks [5S]. 
It can be used to determine the loop-loop scattering am- 
plitudes mentioned in the introduction and visualized in 
Fig.[T] 

For the energy dependence we have introduced in the 
model [31 , 32 , 33 , 16J the two-pomeron scheme of Donnachie 
and Landshoff 15j. Small dipoles couple to the hard and 
large dipoles to the soft pomeron. The transition radius 
was determined through the investigation of the proton 
structure function. Again we refer to the literature for 
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more information and collect some details and the rele- 
vant parameters in the Appendix. 



2.3 Numerical results 

for the nonperturbative corrections 



In order to exhibit the importance of the nonperturbative 
contributions we display their effects explicitly in this sub- 
section. We hasten to add that for the final calculations as 
given in Sec. |3]we do not split our results into perturba- 
tive and nonperturbative parts, but give directly the full 
theoretical results. 



The correction factors 7]^ ( 10 1 and rjx ( 16 1 are displayed in 



Fig. pjas functions of -|-4mj for several vector mesons. 
The lactor 77^ reflects the effect of the distribution of the 
longitudinal momentum; we notice that it remains of or- 
der 1, and its influence is rather weak. In contrast, the 
correction in the transverse amplitude due to the factor 
r]T, reflecting mainly relativistic corrections to the trans- 
verse wave function, is very important, especially for high 
values of + rnj . Its large values at high for the pro- 
duction of p-mesons indicate that the measurement of the 
ratio of the longitudinal to the transverse production cross 
sections tests mainly the wave function. According to our 
calculations, only the quantity rjT for light mesons is very 
sensitive to the specific choice of the wave function. 

There are large effects due to the finite extensions of the 
mesons. We denote by E the ratio of the full amplitudes 
^ to the purely perturbative ones (11 1,(17 1. The ratio E 
is displayed in Fig.|3] left, the results for both polarisations 
are very similar. 

The effects are by no means negligible, even for the heavy 
mesons. The suppression due to the finite extension at low 
values of leads, in the full range of presently available 
data, to a much weaker decrease in than inferred from 
purely perturbative expressions. 

For large values of the quark- antiquark separation Ri, the 
dipole cross section differs from the pure i?^-behaviour as 



given by ([5]),(|12|) 

amplitudes oiT 



We denote by D the ratio of production 
le more realistic dipole cross section ob- 
tained with the stochastic vacuum model divided by the 
result obtained with the purely quadratic expression 

light mesons at low values of Q^. Also D is similar for the 
longitudinal and transverse polarisation. 

The importance of nonperturbative contributions to vec- 
tor meson photoproduction has been stressed before [341 . 
using a different approach to nonperturbative corrections. 



3 Theoretical results 

and comparison with experiment 

A detailed description of the framework of our calcula- 
tions can be found in our previous papers [T6 l [T7 l l2 H [T8| . 



Fig. 2. Correction factors rjL (10 1 and r^T, (16 1 due to the 



longitudinal momentum distribution in the longitudinal and 
transverse wave function of the mesons. The factors are simi- 
lar for the BL and BSW wave functions, except for rjT for light 
mesons. The differences in these cases are illustrated by the 
comparison of the long-dashed and short-dashed curves corre- 
sponding to the BL (long-dashed) and the BSW (short-dashed) 
wave functions for the p-meson. 
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Fig. 3. Correction factors in the evaluation of amplitudes for 
electroproduction. The factor E in the left is due to the finite 
extension of the meson wave function and the correction factor 
D in the right is due to nonperturbative corrections to the 
simple i?i dependence of the dipole cross section. The factors 
are very similar for longitudinal and transverse polarisations. 



The basis for the determination of the loop-loop scatter- 
ing amplitude represented in Fig. [T] is the model of the 
stochastic vacuum. The calculations are determined by 
two parameters, the correlation length of the nonpertur- 
bative gluon fluctuations and the strength of the gluon 
condensate, which can be extracted from lattice calcula- 
tions |5D] . 

All parameters have been determined from other sources, 
and they are collected in Appendix A. Our calculations 
contain therefore no adjustable parameter and all theo- 
retical results are true predictions. 

In the following we present a comparison of our theoret- 
ical results with the available data of elastic photo and 
electroproduction of the S-wave vector mesons, namely 



3.1 Photo and electroproduction of the J/?/; meson 

In Fig. |4] we show the data for the integrated elastic cross 
section a of J/'0 electroproduction at the fixed energy 
W=90 GeV as a function of the photon virtuality . The 
data are from Zeus [551155] and HI [37] collaborations. The 
theoretical calculations (solid line) are made using the BL 



(Brodsky-Lepage) form of the vector meson wave function; 
the BSW-wave function yields very similar results. The 
agreement between data and the theoretical calculations 
is remarkable. The dotted line in the figure, which nearly 
coincides with our theoretical calculation is a fit to the 
data proposed in the experimental paper, with the usual 
form 

" = (T+OVM^ ■ ^^^^ 



Y* p — > J/ij/ p 




W=90 GeV 

- HI 2005 
• Zeus 2004 
o Zeus 2002 



solid: theoretical calculation 

dotted: fitting by experimentalists 
nb 



(GeV^) 



pdf 

Fig. 4. Integrated elastic cross section of J/?/) photo and elec- 
troproduction at the energy W=90 GeV as a function of Q'^. 
The data are from the Zeus |35ll36j and HI [37j collaborations. 
The solid line represents our theoretical calculation using the 
BL wave function. The dotted curve is a pure fit to the experi- 



mental points of the form given by Eq.( 18 1 



Fig. [5] shows the ratio R of the longitudinal to the trans- 
verse cross sections 



R 



(19) 



again for J/t/i elastic electroproduction at W=90 GeV. 
As the ratio cancels influences of the speciflc dynamical 
model, this quantity tests directly details of the overlaps 
of wave functions, helping the study of their longitudinal 
and transverse structures The data are from Zeus [381136] 
and HI 39j,37^ collaborations at HERA. As seen in the 
figure, the presently available data are not very accurate. 
The theoretical calculations with the two kinds of wave 
function - BL and BSW - give nearly the same results for 



(7 = (7^ + (T"^ and for R - 



Fig. [6] shows the comparison of our calculations for the 
energy dependence of cross sections with the recent Zeus 
and HI data [551I351I37] and older photoproduction data 
from fixed target experiments [40(141] at lower energies. 
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Fig. 5. Ratio R of longitudinal and transverse cross sections 
for J/t/j electroproduction as function of for fixed energy 
W=90 GeV. Data from Zeus .38,36 and HI ,39,37^ collabo- 
rations at HERA. The solid line and dashed lines show our 
theoretical calculations respectively with the BL (solid) and 
the BSW (dashed) wave functions. A good numerical repre- 
sentation for the BL result is i?(Q^) ~ 0.76 ((5^/M^)/(l -f 
QVm2)0-09. 



For illustration of our results for electroproduction we 
show data and theoretical results in the range 6.8 - 
7 GeV^ for which there are both Zeus and HI data. 

Often the W dependence of experimental cross sections is 
fitted through single powers in the form 



Const. X M^-'C^") 



(20) 



which may be useful in limited W ranges. Values of the 
parameter 5 obtained at several values of in J/?/; elec- 
troproduction are compared to our theoretical predictions 
in Fig.|7] based on the two-pomeron scheme. 

In contrast to purely perturbative approaches, our theo- 
retical treatment allows to calculate the dependence of the 
cross section on the momentum transfer t. The data on the 
^-distribution in J/V' photo |35|l37j and electroproduction 
[551157] at W=90 GeV are shown in Fig. [s] together with 
our theoretical results (solid lines). Also shown are elec- 
troproduction data and calculations at — 6.8 — 7 GeV^ 
where measurements from both Zeus and HI are available. 

Our calculation predicts a curvature in the log plot for the 
t distribution, which we may describe with the form 



da 

For 7*p - 
give a 



da 



d\t\ 



xF{\t\)^ 



da 
d\t\ 



-m 



J t=0 



{i + a\t\y 



(21) 



pip at W=90 GeV our calculations at — 
■ 4.06 GeV~2 and b = 1.75 GeV"^ Present 



_ y* p ^ J/f p / 


Q^=6.8-7.0 GeV^ 














- HI (2005) 


ZEUS (2004) 



100 1 50 

W (CieV) 



Fig. 6. Energy dependence of the integrated elastic cross 
section for J/tp production from the HERA collaborations |35l 
36,37 , compared with our theoretical calculations. The fixed 
target photoproduction data are from the E-401 and E-516 
experiments |40II41| . Electroproduction is represented by the 
U range 6.8-7.0 GeV^ for which there are both Zeus and HI 
data points. 



data give no clear cut evidence for this curvature. The 
distribution becomes flatter as increases, but in the 
investigated range does practically not change with the 
energy (no shrinking). 

In order to determine the integrated cross section the ex- 
perimental data have to be extrapolated to the point of 
minimum momentum transfer. This brings a theoretical 
bias to the experimental points. It is therefore meaningful 
to compare our theoretical results with observed data at 
small but finite momentum transfer. This is done in Fig. 
|9] the agreement is excellent. 
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Fig. 7. Q dependence of the 5 parameter describing the 
energy dependence of the J/t/j integrated elastic cross section, 
from Zeus and HI collaborations |35ii36ii37, . compared to our 
theoretical predictions. 



3.2 Photo and electroproduction of T meson 

For the T there are only two data points for photoproduc- 
tion, at <W >=120 GeV ggjand < W >=143 GeV gg, 
shown together with our theoretical energy dependence in 
In the second plot we show our calculation of the 
dependence at the fixed energy W=130 GeV, together 



Fig 
n2 



with the data for = at the energies 120 and 143 GeV. 
Both plots show a fairly good agreement with the data. It 
should be noted that, in spite of the high mass scale, the 
finite extension of the meson has a considerable effect, as 
can be seen in Fig.[3]left. 



3.3 Photo and electroproduction of p meson 




7* p ^ J/V|/ p 

W=90 GeV 

Q2=6.8 GeV^ 



- Zeus 2004 

- HI 2005 

solid: theoretical calculation 
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Fig. 8. The t dependence of the differential cross sections of 
J/t/i elastic photo 35,37 and electroproduction 36i32]- The 
Zeus(2002) photoproduction results have separate information 
from events with decaying into fi'^ j-L~ (full circles) and e'^e" 
(empty circles). The solid lines represent our theoretical re- 
sults. 



The classical fixed target experiments |H1H51H6II47IH5] pro- 
vide important reference for the magnitudes of cross sec- 
tions in p photo and electroproduction at center of mass 



energies about 20 GeV. The data are shown in Fig 
together with the results of our calculation. The agree- 
ment is satisfactory given the quite important discrepan- 
cies within the data in the interval — A — 7GeV^. 

Fig. [12] shows the data for the integrated elastic cross sec- 
tion a of p electroproduction at the fixed energy W=90 
GeV as a function of the photon virtuality . The data 
are from Zeus [38,49 and HI [SUIIST] collaborations. The 
solid line represents the theoretical calculations using the 
BL form of the vector meson wave function. Our theoreti- 
cal calculations give an good overall description of the data 
but do not reproduce the details of the dependence of 



p electroproduction as well as they did in J /ip production. 
The present data indicate that in the p meson case a fit- 
ting of a simple form like Eq. (18) is not satisfactory for 



11 the whole range. Two of such forms may be applied 



with a transition in the values of both parameters (normal- 
ization and power) occurring somewhere in the range 
from 5 to 7. More experimental measurements are needed 
to clarify the structure of the data in this region, which 
may contain important information about the dynamics 
of the process. 



Fig. 
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shows the ratio R = jG^ of cross sections for 
longitudinal and transverse polarisations, for p elastic elec- 
troproduction at W=20 and 90 GeV. The data at 90 GeV 
are from NMC gS], Zeus [351113 and HI [SHIMJISI]- The 
W=20 GeV data are from the E-665 experiment [47|. To 
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Fig. 9. 



and W dependences of the (nearly) forward dif- 



ferential cross sections of J/i/) elastic photo 35 ,37 and elec- 
troproduction |36II37) . The solid lines represent our theoretical 
results. 



indicate the differences, the theoretical results for both 
BL (solid line) and BSW (dashed line) wave functions are 
displayed. For large Q^, R becomes very sensitive to the 
values of the small transverse cross section. The measure- 
ment of this quantity R therefore provides important tests 
on the structure of the meson wave function. 



Fig. 14 shows the energy dependence of p electroproduc- 
tion for several values of together with our theoretical 
results. On the right-hand-side plot we show the effective 
power 5 for the energy dependence (see eq. (20)) for the 
M^-region of approximately 30 to 130 GeV. The experi- 
mental points are from Zeus 38,49,53] and HI |B0,51 . 
Our theoretical description, solid line, is very satisfactory 
within the experimental errors. 

There are no published measurements of the differential 
cross section da/d\t\ in p production for nonzero values of 
. The Zeus photoproduction data at W=75 and 94 GeV 



Fig. 10. Integrated elastic cross section of T photoproduction 
as function of the energy and electroproduction at the energy 
W=130 GeV as a function of . The data at the energies 
< VK > = 120 and 143 GeV are respectively from Zeus 03 s^nd 
HI [43| collaborations. The solid lines represent our theoretical 
calculations using the BL wave function. 



are shown in Figs. 15 In the low t range there are 
new preliminary data from HI [55] at several energies. The 
numbers for the differential cross sections at 70 GeV for 
low t extracted from their plots are included in the figure, 
and they seem to confirm the previous Zeus measurements 
at W=75 GeV [53]. Our theoretical calculations are shown 
in the figures, the description is satisfactory. 



3.4 Photo and electroproduction of u meson 



As particles of about the same size and mass, u and p 
have similar behaviour in the soft processes that we study 
here. 
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Fig. 11. Integrated elastic cross section for p electroproduc- 
tion at the energy W=20 GeV as a function of the photon 
virtuaUty . The data are from fixed target experiments [441 
I45ll46ll47ll48j . spanning almost two decades. The solid line rep- 
resents the results of our calculations. 




Fig. 12. Integrated elastic cross section of p electroproduction 
at the energy W=90 GeV as a function of Q^. The data are 
from Zeus |38ll49j and HI 50,51 collaborations. The solid line 
shows our theoretical results using the BL wave function. 




(GeV') 



Fig. 13. Ratio R between longitudinal and transverse cross 
sections for p electroproduction as function of for fixed 
energies W=20 GeV , with data from E-665 gT], and W=90 
GeV with data from NMC [48], Zeus [38ll49] and HI [5211501151] . 
The solid and dashed lines show the theoretical calculations 
with the BL and BSW wave functions respectively. 

The existing data [551155] on the t dependence of the dif- 
ferential cross section in lo production are shown in Fig. 
pT] In the figure are put together the data points of the 
energies W=15 GeV (with ^ q) and W=80 GeV (with 
Q^=0.1 GeV^), and the corresponding theoretical curves. 
The coincidence of shapes exhibits the universality of the 
form factors of t dependence in our model. Our prediction 
of a curvature in the plot of da j d\t\ seems to be confirmed 
by the data. 



The dependence of w photo and electroproduction f56', 
157] is shown in Fig. 16 The agreement between our results 
and the data is not perfect, but satisfactory, in view that 
there is no free parameter involved in the calculations. 



3.5 Photo and electroproduction of meson 

The (/) meson has a strategic place between the J/?/; and 
the p meson and may help to understand the differences 
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Fig. 14. Energy dependence of p electroproduction cross sec- 
tion in the region from approximately 30 to 130 GeV. The 
solid lines show our theoretical calculations for several values 
of and the data are from Zeus [38ll49ll53| and HI [SOllST] 
collaborations. The plot in the right hand side shows data and 
theoretical values for the parameter 5 of the energy dependence 



Fig. 15. Experimental data and our theoretical prediction 
(solid line) for the t dependence of p photoproduction cross 
sections. The published data are from the Zeus collaboration, 
at the energies 75 GeV (low t) [53] and 94 GeV (large t) [54]. 
The 94 GeV data are rescaled with a factor (90/94)°-^® = 0.993 
in the figure. We also include in the low-t figure preliminary 
information from Ifl [55j . extracted from their plots, for W=70 
GeV. 



in behaviour of heavy and Hght vector mesons and also to 
clarify the interplay of perturbative and nonperturbative 
aspects of QCD. 



In Fig. 



18 we show the dependence of the integrated 
photo and electroproduction at 



elastic cross section in 
W=75 GeV [60.i61ll62] . together with our theoretical cal- 
culations. As in the case of p electroproduction, there is 
indication that the data cannot be w ell represented by 
a single expression of the form of Eq. (18 1. The data are 



poorer here than in the p case, and it is important to inves- 
tigate the possibility of a transition region in an interme- 
diate range below 10 GeV^ in which the normalization 
and the power change values rather rapidly. 



The data for the ratio R = ja^ for W=75 GeV as func- 
tion of [6l][62| are shown in Fig. 19 together with our 
results with BL and BSW wave functions. The calcula- 
tions exhibit, as in the p meson case, the sensitivity of the 
ratio R to details of the wave functions. Offering a refer- 
ence for the two kinds of calculation, the plot also shows 
(dotted line) a fit of the data, made by experimentalists. 

The effective power 6{Q^) describing the energy depen- 
dence as in electroproduction [62j is shown in Fig. 
|20| and compared with our results using the two-pomeron 
model. The theoretical dependence is similar to that of 
J/iJj and p production, whereas the experimental points. 
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Fig. 16. dependence of the cross section of uj elastic pro- 
duction |56II57| . Tlie solid line represents our calculations. 



Fig. 18. dependence of production cross sections at 

W=75 GeV and our theoretical description. The data are from 
Zeus |60ii62) and HI [51]. 




Fig. 17. t dependence of u photoproduction cross sections [581 
159] . The plot puts together data points at the energies W=15 
GeV (with g2 = 0) and W=80 GeV (with g^^O.l GeV^), 
and the corresponding theoretical curves. The curvature and 
the similarities of shapes of the t distributions for two differ- 
ent energies and values are characteristic features of our 
framework. 



Fig. 19. dependence of the ratio R = ja^ in electro- 
production [611162] at fixed energy W=75 GeV. The solid and 
dashed lines represent our calculations with BL (solid) and 
BSW (dashed) wave functions; the dotted line is a fit of the 
form R = 0.51 (QV^^)"'**^- 



with large errors, indicate a flatter behaviour. Obviously 
more data are necessary. 

Our calculations give very good descriptions for the 
dependence of the differential elastic cross section in for- 
ward directions in electroproduction, as it does in the calculations give a very satisfactory description. 



J/ip case. Fig. 21 shows the comparison with the experi- 
mental data |62j. 

There are no published measurements of da/d\t\ in pro- 
duction for nonzero . The Zeus photoproduction data 
at 94 GeV [MlISi] are shown in Fig. [22] Our theoretical 
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Fig. 20. The effective power 5 parameter governing the energy 
dependence of (j) electroproduction. The data are from Zeus 
[52] . the sohd line is our theoretical result. 



3.6 Results concerning several vector mesons 

The quantitative predictions made in a unique way for dif- 
ferent kinds of vector mesons cover several scales of mag- 
nitudes in cross sections. This global coverage is exhibited 
in Fig. 23 The same global description is given for the 
forward differential cross section shown in Figj24] In this 
figure the charge factors squared e^, see Eq. (jsk, for each 
kind of meson are extracted, making the quantities almost 
universal in a -|- My plot. However, the universality is 



only approximate, and our calculation predicts correctly 
the observed displacements. 

The data on the energy dependence of the integrated cross 
sections for p, ^ and mesons have been presented and 
compared to the theoretical predictions for each case. The 
parameter 5{Q'^) of the suggested simple energy depen- 
dence 

CT(g2) ^ Const. X 1^5(0') 

has also been given in each case. This parametrization is 
an approximation valid in a limited energy range, since the 
true energy dependence in our calculation is determined 
by the two-pomeron scheme, but it is considered useful in 
practice. 

We then evaluate 5 using the energy range W=20 - 100 
GeV, for all values of . The results are put together 
in Fig. |25| We note that all curves start at the minimum 
value 4 x 0.08 at the same unphysical point + My = 
and all are asymptotic to (5 = 4 x 0.42 as increases. We 
then have the form of parametrization 



5{Q^) = 0.32 -f- 1.36 



(1 



'/Ml 



A + {1 



(22) 




Fig. 21. dependence of the differential cross sections of 
elastic in the forward direction. The solid lines represent our 
calculation. The second plot contains the same information, 
showing the typical behaviour of a straight line in the variable 



Meson 


P 







r 


A 


124.926 


51.9183 


2.0624 


0.9307 


n 


1.239 


1.239 


1.12 


0.22 



Table 1. Values of A and n in the exponent S{Q^) from Eq. 
(|22l) 



and the values for A and n are given in table (T) 

Our model has definite predictions for the t dependence in 
differential cross sections. The shape can be conveniently 
represented by the form given in Eq.(21l. Fig. 26 shows 
the form factor for all vector mesons at fixea energy 

W=90 GeV. 

The curvature in the log graph is a prediction of our frame- 
work. The values of the parameters as functions of are 
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Fig. 22. t dependence of photoproduction cross sections and 
its theoretical description in our nonperturbative calculation 
with the stochastic vacuum model. The published data are 
from the Zeus collaboration, at the energy 94 GeV, first at low 
t |60| and then at larger t [54j . 



Fig. 24. Differential cross sections in a forward direction for 
all vector mesons at W=90 GeV, as functions of +My, with 
extraction of charge factors . The lines represent the theoretical 
calculations with the stochastic vacuum model as in Fig. |23[ 




Fig. 23. Integrated elastic cross sections for all vector mesons 
at W=90 GeV, as functions of + My . The lines represent 
our theoretical calculations. 



shown in one of the plots. In Table |2] some numerical val- 
ues of a, 6 are given. 

As the virtuality grows, the ranges of the overlap func- 
tions decrease, and the electroproduction cross sections 
of all mesons become flatter, all tending together to the 
shape characteristic of the T meson, with same limiting 



Fig. 25. Parameter S{Q'^) of the energy dependence of the 
cross sections. 



values a — 4.02 GeV^ and b — 1.60 GeV^ for the parame- 
ters. The limiting shap e of the distribution for very large 
is illustrated in Fig. 27 where we see all vector mesons 
superposed. In the figure we draw the bit of straigth line 
representing the slope considered as the average for the 
interval from \t\ = to \t\ = 0.2 GeV^. As indicated in- 
side the plots, the calculations of form factors presented in 
the figures are made for W — %Q GeV. In the second plot 
presented in Fig. 27 we show the (absence of) dependence 
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Fig. 26. Form factor -F(|tj) of the \t\ distribution in the 
elastic differential cross sections photoproduction of vector 
mesons, with the characteristic curvatures in the log plot pre- 
dicted i n o ur calculations. The curves follow the shapes given 
by Eq.(21l, with a{Q^) and b{Q^) represented in the sec- 
ond figure. Specific numerical values of the parameters for 
= , 10 and 20 GeV^ are given in Table [2] The val- 
ues in the limit of very large are a = 4.02 GeV~^ and 
b = 1.60 GeV-^ 



of the form factor on the energy W. Thus we predict that 
there is no shrinking of the forward peak in da/d\t\ as 
the energy increases. The experimental data are not yet 
sufficient to test all these predictions. 



4 Relation with other approaches 

Our work runs parallel to other developments , of mainly 
nonperturbative nature, based on the colour dipole dy- 
namics. These treatments, which started with studies of 
7*p total cross section [63] and soon were extended to 



Table 2. Some values (in GeV ^ ) of the parameters a and b of 
X 



dtr 
d\t\ 



t=0 

is shown in Fig. [26| 



(l + a|t|)^ 



in Eq.(21 \. The full Q dependence 
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10 GeV^ 




20 GeV^ 


Meson 
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a 
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p(770) 


7.06 


3.09 


4.349 


1.996 


4.145 


1.846 


a;(782) 


7.20 


3.08 


4.323 


1.990 


4.141 


1.835 


0(1020) 


5.40 


2.77 


4.211 


1.934 


4.091 


1.807 


j/V(is) 


4.06 


1.75 


4.026 


1.696 


4.022 


1.670 


r{is) 


4.03 


1.61 


4.027 


1.606 


4.025 


1.604 



diffractive deep inelastic scattering [M', are based on the 
assumption that the photon fluctuates into a qq pair, which 
then scatters with the proton, according to a properly 
built dipole cross section CTdipoic- The qq dipole then re- 
combines to form again a photon, or eventually a vector 
meson, in the final state. Photons, either real or virtual, 
enter in a well defined way through their QED wave func- 
tions, and the dipole cross section is built with the neces- 
sary ingredients to describe the observed phenomenology 
and, as much as posible, to follow QCD prescriptions. The 
proton structure is introduced to give an impact parame- 
ter dependence that allows description of t dependence in 
a predictable way. 

Concerns with observed saturation effects in the and 
X dependences led to an ansatz for (Tdipoic of the form 



CTdipoic (a;, r) = o-o(l - exp (-r^Q^(a;)/4)) , 



(23) 



where ctq is a constant and Qs{x) denotes x dependent 
saturation scale. For processes governed by quarks of mass 
m/ , the Bjorken scale x is taken as 



(24) 



This kind of dipole cross section worked well for what it 
was aimed for, giving a description of the HERA data on 
inclusive (related to F2) and diffractive 7*p processes. In 
an effort to improve the model including scaling violations 
and QCD evolution, B artels, Golec-Biernat and Kowalski 
[55] replaced the saturation scale by a gluon density 
obeying DGLAP evolution. The interplay between satura- 
tion and evolution improved fits of the F2 data , especially 
for large Q^. 

Another improvement of the original saturation model 
came with the Color Glass Condensate model [66 , dealing 
with the question of the low- a; behaviour near the satura- 
tion region through BFKL dynamics. These calculations 
were successful as an extension of the model, although 
they did not include charm contributions. 

An important line of development of the phenomenolog- 
ical color dipole models came with the introduction of 
effects of the proton shape, with explicit impact parame- 
ter dependence in the dipole scattering, by Kowalski and 
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Fig. 27. The plot in tiie left hand side shows the shape of the 
|t| distribution in the differential cross section for electropro- 
duction common to all vector mesons for very large . The 
straight line passes through the points \t\ equal to zero and 0.2 
GeV^, indicating the average slope that would be measured 
in this limit. The plot in the right hand side shows that J'(ltl) 
does not depend significantly on the energy, and thus measure- 
ment of slopes at fixed would give a constant value for each 
vector meson. 



Teaney ^7j, and consequently with possibility of treat- 
ment of t distributions. 

This work was extended by Kowalski, Motyka and Watt to 
the description of photo- and electroproduction of vector 
mesons going beyond F2 and DVCS. The authors 
use the same kind of dipole cross section as KT , based 
on Glauber-Mueller: 



lution without quarks. Application of this model requires 
a number of ingredients and parameters, as form of gluon 
density at ^q, effective x for heavy quarks, proton shape 
T(&), b-dependence in the form of Udipoie- The paper de- 
scribes part of the recent HERA data on photo- and elec- 
troproduction of vector mesons that we discuss in the 
present paper, with interesting similarity of results. 

A fundamental and original approach for the implemen- 
tation of the dipole treatment of vector meson produc- 
tion, which provides the framework for the present work, 
was developed by the Heidelberg group pT |[24ll25] . This 
approach calculates the basic loop-loop interaction using 
the stochastic vacuum model [T31[TH[5D], which is a gen- 
uine nonperturbative treatment that allows connections 
with other branches of nonperturbative QCD, especially 
with the fundamental and striking feature of confinement. 
The method is based on functional methods [TT], and has 
been used recently for investigations on the foundations 
and limitations of the dipole model [T2l[69ll70] . 



The same approach using the loop-loop ampUtude and the 
stochastic vacuum model was used consistently by Don- 
nachie and Dosch for calculations of DVCS [55] and struc- 
ture functions [33 . 

The treatment of photoproduction of all vector mesons 
through a dipole Pomeron model respecting the unitarity 
bounds by Martynov, Predazzi and Prokudin |71j nicely 
described W and t dependences in elastic photoproduc- 
tion processes, and was soon extended to electroproduc- 
tion 

A separation of hard (small) and soft (large) dipole inter- 
actions, testing several prescriptions for the dipole cross 
section and several forms of wave function, was used by 
Forshaw, Sandapen and Shaw [73|[74] to describe energy 
and dependences of the total j*p cross section, diffrac- 
tive electroproduction, DVCS and exclusive J/ip electro- 
production. The separation of dipoles in two classes is 
similar to the two-pomeron model adopted in the present 
work. 

As a whole, these nonperturbative models for 7* induced 
processes, all based in the dipole picture, are phenomeno- 
logically satisfactory. This is not surprising, since it has 
already been shown [T5] that many features of these pro- 
cesses, particularly the dependence, are reproduced by 
the overlap of the light cone wave functions of photons and 
mesons folded with the basic behaviour of the dipole 
cross sections. 

Unfortunately the presently available data are not suffi- 
cient to discriminate details of different approaches, par- 
ticularly in the incorporation of energy dependence. 



da, 



dipole 



2[1 



exp 



27V, 



■r^as{^l'^)xg{x,^J.'^)T{b) 



(25) 

The scale jj? is related to the dipole size, and the gluon 
density is evolved from /1q to using LO DGLAP evo- 



5 Summary and discussion 



We have shown the predictions for elastic electroproduc- 
tion processes using two basic ingredients: 
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1) the overlaps of photon and meson wave functions, writ- 
ten as packets of quark-antiquark dipoles, with protons 
described also as packets of dipoles (in a convenient di- 
quark model for the nucleon), 

2) the interaction of two dipoles described in terms of ge- 
ometric variables in an impact parameter representation 
of the amplitudes based on nonperturbative properties of 
the QCD gluon field. 

These quantities put together and integrated over the dis- 
tribution of dipoles in initial and final states lead to a cor- 
rect description of the data concerning all vector mesons. 

In all cases the variations with energy are very well de- 
scribed by the Regge picture, with soft and hard pomerons 
coupled to large and small dipoles, respectively. We recall 
that in approaches mainly based on perturbation theory 
the energy (or x) dependence is introduced through the 
gluon distribution in the proton. 

Each of the different mesons enter in the calculation char- 
acterized only by the masses and charges of its quark con- 
tents, and with their normalized wave function individual- 
ized only by the corresponding electromagnetic decay rate 
(related to the value of the wave function at the origin) . 

The specific nonperturbative input is the stochastic vac- 
uum model, which has been successfully applied in many 
fields, from hadron spectroscopy to high energy scatter- 
ing. The basic interaction of two dipoles depends only on 
universal features of the QCD field, which are the numer- 
ical values of the gluon condensate and of the correlation 
length of the finite range correlations. These two quan- 
tities have been determined by lattice investigations and 
tested independently in several instances of phenomeno- 
logical use of the dipole-dipole interaction. 

As has already been pointed out [TS] , the main features of 
the dependence of electroproduction of vector mesons 
are contained in the overlap integral. For values of at- 
tainable at present this overlap integral is determined by 
perturbative QCD, through the photon wave function, and 
by nonperturbative QCD, through the meson wave func- 
tion. The production of transversely polarized mesons is 
determined by the meson wave function even for very high 
values of Q^. The importance of nonperturbative effects 
in vector meson production has also been stressed in [51] . 
where the perturbative two gluon exchange has been sup- 
plemented by an exchange of nonperturbative gluons. In 
our approach, which incorporates both perturbative and 
nonperturbative effects, we are able to give a fair overall 
description of all observables of vector meson production: 
the energy dependence, the dependence, the ratio of 
longitudinal to transverse mesons and the angular distri- 
bution. 

In this paper our calculations are compared to the large 
amount of HERA data for p,uj, 0, J/ip and T mesons. 

Botli autliors wish to thank DAAD (Germany), CNPq (Brazil) 
and FAPERJ (Brazil) for support of the scientific collaboration 
program between Heidelberg and Rio de Janeiro groups work- 



ing on hadronic physics. One of the authors (EF) is grateful to 
CNPq (Brazil) for research fellowship and grant. 



A QCD and wave function parameters 

Here we recall some expressions and the numerical values 
of quantities used in the nonperturbative contributions, 
which have been unchanged for many applications. 

The loop- loop scattering amphtude S{b,W, zi, Ri, Z2 — 
1 /2, i?2) determining the essential quantity J{q, W, zi,Ri) 
in Eq. (jsj) , can be calculated using the stochastic vacuum 
model. The essential input parameters of this model are 
the correlation length a of the gauge invariant two gluon 
correlator and the gluon condensate (g'^FF). These quan- 
tities have been determined in lattice calculations [5U] . 
The numerical values used in this and previous papers are 



= 0.346 fm 



(g^i^F) a'^ = 23.5 



(26) 



The constant k — 0.74 appearing in the correlation func- 
tions preserves its value determined by lattice calculations. 

The energy dependence is based on a two-pomeron model, 
small dipoles with size R < Rc couple to the hard pomeron 
with an intercept aph{0) — 1.42, whereas large dipoles 
with R > Rc couple to the soft pomeron with an intercept 
a_ps(0) = 1.08. The transition radius Rc has been de- 
termined from the x-dependence of the proton structure 
function to i?c = 0.22 fm. 

For the proton wave function occurring in Eq. (§ we use a 
diquark-quark Gaussian wave function with the transverse 
radius Rp = 0.75 fm. 

The light quark masses which can simulate confinement 
effects in the otherwise perturbative photon wave function 
[2T are 



niu — md — 0.2 GeV 



= 0.3 GeV 



(27) 



For the heavy quarks we take the renormalized masses in 
the MS scheme 



= 1.25 GeV 



nib = 4.2 GeV (28) 



The full form of photon and vector meson wave functions 
used in our work, including the helicity dependences, have 
been presented before [W.fTB], with two forms for the vec- 
tor mesons: the Bauer-Stech-Wirbel (BSW) [26 , and the 
Brodsky and Lepage (BL) j27j wave functions, where the 
separate r, z dependences are respectively 




(29) 



zz exp 



2a;2 



exp[ 



(here My represents the vector meson mass) and 




(30) 



exp 



1N2 
2' 



IzzuP- 



exp[ 



-2zzuP'r'^ 
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Table 3. Parameters of the vector meson wave functions 









BSW 






BL 










transverse 


longitudinal 


transverse 


longitudinal 


Meson 


fv (GeV) 


a;(GeV) 


N 


tj(GeV) 


N 


c^(GeV) 


N 


c^(GeV) 


iV 


p(770) 


0.15346 


0.2159 


5.2082 


0.3318 


4.4794 


0.2778 


2.0766 


0.3434 


1.8399 


(j(782) 


0.04588 


0.2084 


5.1770 


0.3033 


4.5451 


0.2618 


2.0469 


0.3088 


1.8605 


0(1020) 


0.07592 


0.2568 


4.6315 


0.3549 


4.6153 


0.3113 


1.9189 


0.3642 


1.9201 


J/^(1S) 


0.26714 


0.5770 


3.1574 


0.6759 


5.1395 


0.6299 


1.4599 


0.6980 


2.3002 


T{1S) 


0.23607 


1.2850 


2.4821 


1.3582 


5.9416 


1.3250 


1.1781 


1.3742 


2.7779 



with ruf representing the quark mass. 

The parameters N (normahzation) and lu (that fixes the 
extension) are determined using the electromagnetic decay 
rates of the vector mesons. Their values are collected [TB] 
in Table M 

After summation over helicity indices, the overlaps of the 
photon and vector meson wave functions 

P7*v,a((9^;^;i,R.i) = il:vx{zi,'R.i)*'ip'y'x{Q^; zi,'Ri) (31) 

that appear in Eq.([2]) are given by 



P'y'V,±i-BSwiQ ;z,r) = ev- 



2n 



(32) 



(^ef uj'^r[z^ + z^]Ki{ef r) + mjKo{ef r)j 4>BSw{z,r) 



and 



P'y'V.,±i:BL{Q'^;z,r) = ev 



V6a 



(33) 



for the transverse case, BSW and BL wave functions re- 
spectively. For the longitudinal case we can write jointly 



P-i'V,0;x{Q ]z,r) ^ -IQev— — c^z z Q Ka{ef r)(t)x{z,r) , 

(34) 

with X standing for BSW or BL. The effective quark 
charges ey are l/\/2 for p, l/3\/2 for -1/3 for (j) and 
T and 2/3 for ill. The quantity 



^z{l - z)Q^ + mj 



(35) 



and the modified Bessel functions are introduced by the 
photon wave functions. 
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